Ab initio calculations of third-order elastic constants and related properties for selected 

semiconductors 
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We present theoretical studies for the third-order elastic constants Cijt in zinc-blende nitrides A1N, 
GaN, and InN. Our predictions for these compounds are based on detailed ab initio calculations of 
strain-energy and strain-stress relations in the framework of the density functional theory. To judge 
the computational accuracy, we compare the ab initio calculated results for djk with experimental 
data available for Si and GaAs. We also underline the relation of the third-order elastic constants 
to other quantities characterizing anharmonic behaviour of materials, such as pressure derivatives 
of the second-order elastic constants cL and the mode Griineisen constants for long-wavelength 
acoustic modes 7(q,j). 
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I. INTRODUCTION 

Third-order elastic constants are important quan- 
tities characterizing nonlinear elastic properties of mate- 
rials and the interest in them dates back to the beginning 
of modern solid state physics j 1 ' 2 ' 3 ' 4 ' 5 ' Third- and higher- 
order elastic constants are useful not only in describing 
mechanical phenomena when large stresses and strains 
are involved (e.g., in heterostructures of optoelectronic 
devices), but they can also serve as a basis for discus- 
sion of other anharmonic properties. The applications 
include phenomena such as thermal expansion, temper- 
ature dependence of elastic properties, phonon-phonon 
interactions etciS 

As far as theoretical studies are concerned, at the be- 
ginning the third-order elastic constants were calculated 
in the framework of the valence force Keating model— 
Later on, many other more sophisticated microscopic 
theories were employed to describe and predict nonlin- 
ear elastic properties of crystals on the basis of their 
atomic composition.^ Nowadays, precise ab initio calcu- 
lations seem to be the most promising approach to handle 
this task. Such applications of density functional theory 
(DFT) on the local density approximation level (LDA) 
were already reported^ 9 - 

Recently, one observes increased interest in nonlinear 
effects in elasti c 10 : 11 ' 12 and piezoelectric properties! 13 : 14 
This is strongly connected to the fact that research fo- 
cuses nowadays on the semiconductor nanostructures. 
In such systems these nonlinear effects are not only 
more pronounced than in bulk materials, but very of- 
ten their reliable quantitative description is a prerequi- 
site for correct theoretical explanation of the experimen- 
tal data j 12 : 14 ' 15 : 16 ' 17 : 18 In this paper, we perform ab initio 



calculations of the unknown third-order elastic constants 
in cubic nitrides. The nitrides are technologically im- 
portant group of materials for which the nonlinear ef- 
fects are particularly significant i 10 : 12 ' 19 : 20 Therefore, the 
knowledge of the third-order elastic moduli will definitely 
improve the modeling of nitride based nanostructures. 
In this work we also briefly discuss the applications of 
Cijk to determination of other anharmonic properties, 
namely, pressure derivatives of second-order elastic mod- 
uli and mode Griineisen constants 7(q,j). Since the 
third-order effects are rather subtle, their computational 
determination can also serve as a precise test of accuracy 
for modern ab initio codes based on DFT approach. 

The paper is organized as follows. In Sec. [IT] we give a 
general overview of the nonlinear elasticity theory. Sec. 
IIIII contains a description of employed methodology. Also 
results for third-order elastic constants obtained from ab 
initio calculations are presented there. Our findings for Si 
and GaAs are compared with previous numerical calcula- 
tions and measurements, later on theoretical predictions 
for zinc-blende nitrides A1N, GaN, and InN are given. 
Sees. IIVI and fVl deal with the determination of quantities 
related to third-order elastic constants, namely, the pres- 
sure dependent elastic constants and mode Griineisen 
constants, respectively. Finally, we conclude the paper 
in Sec. IYT1 



II. 



OVERVIEW OF NONLINEAR ELASTICITY 
THEORY 



Here we will recall some basic facts from nonlinear the- 
ory of elasticityfii 2 -' 3 -' 4 -^ Let us consider point a which, 
after applying strain to a crystal, moves to the position 
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x. After introducing the Jacobian matrix J 

_ dxj 

ddj 

we may define the Lagrangian strain 
1 



(i) 



(2) 



ijkViVjVk 



which is a convenient measure of deformation for an elas- 
tic body. 

The energy per unit mass E(rf) corresponding to the 
applied strain may be developed in power series with re- 
spect to r). This leads to the expression 

p E(rj) = ^ ^ djViVj + 

i,j = lfi i,j,k=l,6 

(3) 

where we applied Voigt convention (7711 — > 771, 7722 — » 
772, %3 -* V3, V23 W 2 , 7713 T75/2, 7712 jye/2) 
and introduced the density of unstrained crystal pq. The 
Cjj and Cijk denote here second- and third-order elastic 
constants respectively^ If we introduce J = (1 + e) and 
assume that e is symmetric (rotation free) linear strain 
tensor, the definition of 77 [Eq. @] yields 

1 



III. DETERMINATION OF THIRD-ORDER 
ELASTIC CONSTANTS 

A. Methodology and computational details 

In this work, we have determined third-order elastic 
constants for Si, GaAs, and zinc-blende nitrides (A1N, 
GaN, and InN) on the basis of quantum DFT calculations 
for deformed crystals. The results were obtained in two 
ways - employing strain-energy formula [Eq. (f5])] and 
from strain-stress relation [Eqs. © and 0]. 

The detailed procedure was as follows. We considered 
six sets of deformations parametrized by 77 



« = e+-e . 



(4) 



Substituting the above result to the expansion in Eq. ([3]) 
and leaving only terms up to second order with respect 
to components of e recover the infinitesimal theory of 
elasticity. 

Naturally, the general expression for energy of strained 
crystal, as given by Eq. ([3]), can be simplified by em- 
ploying symmetry considerations. For cubic crystals, this 
procedure yields the following formula: 

p E(f]) = icn (771 + 772 + -qt) + ic 44 (77I +77e) + 



(5) 



+012(771772 + 773772 + 771773) + 
gCm (77? + v 3 2 + 773) + 

2 C ii2 {v2Vi + V3Vi + V2V1 + V3V1 + V2V3 + V2V3) + 

C123 Vi V2V3 + ^144 (77177! + 77277! + 77377!) + 

^Ci55 {V2V4 + V3V4 + V1V5 + V3V5 + ViVe + V2V&) + 

C4567?4775?76 + . . . 



Another fundamental quantity in the theory of finite 
deformations is Lagrangian stress 

Pott—, (6) 



U 



dm. 



which can be expressed in terms of linear stress tensor cr 
using the following formula 



t = det{J)J- 1 <r{j T ) 1 



Again, Voigt convention (t xl — > t\, t 2 2 - 
h3 — > h, ti 3 — > t 5 , ti2 is used here. 



U, t 



33 



(7) 



Va 

Vd 
Ve 



(77,0,0,0,0,0), 
(77,77,0,0,0,0), 
(77,0,0,77,0,0), 
(77,0,0,0,77,0), 
(77,77,77,0,0,0), 
(0,0,0,77,77,77). 



(8) 



In every case, 77 was varied between —0.08 and 0.08 with 
step 0.008. For every deformed configuration, the po- 
sitions of atoms were optimized and both energy and 
stress tensors were calculated on the basis of quantum 
DFT formalism. In this way, for each type of distortion, 
dependencies of energy ^(77) and stress tensor t(rf) on 
strain parameter 77 were obtained. The numerical results 
have been in turn compared with the expressions from 
the nonlinear theory of elasticity, which are summarized 
in Tablefl] This allows to extract the values of the second- 
and third-order elastic constants, by performing suitable 
polynomial fits. 

The DFT calculations have been performed using the 
ab initio total energy code VASP developed at the Insti- 
tut fiir Materialphysik of Universitat Wieni2Li22i2£ The 
projector augmented wave (PAW) approach^ has been 
used in its variant available in the VASP packaged For 
the exchange-correlation functional generalized gradient 
approximation (GGA) according to Perdew, Burke and 
Ernzerhof (PBE) 26 i 27 has been applied. For Ga and In, 
semicore 3d and 4d electrons have been explicitly in- 
cluded in the calculations. 

Since the determination of subtle third-order effects 
requires high precision, we have performed careful con- 
vergence tests for parameters governing the accuracy of 
computations. On the basis of our tests we have chosen 



the following energy cutoffs Cutoff 



plnN 
^cutoff 



600 eV, E^i - 
800 eV. For the 



700 eV, andi^ t N off = i& 
Brillouin zone integrals we have followed the Monkhorst- 
Pack scheme^, in Si and GaAs we have used 13 x 13 x 13 
mesh, whereas for A1N, GaN, and InN we have applied 
llxllxll sampling. One example of performed tests for 
GaN is presented in Fig. Q] It illustrates the dependence 
of two sample elastic moduli Cm and C144 on the energy 
cutoff and density of Monkhorst-Pack k-point mesh. For 
the chosen parameters (Setoff = ^00 eV and llxllxll 
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TABLE I: Dependencies of energy and stress on deformation parameter r\ for considered types of deformation t]a, ■ ■ ■ ,riF, which 
have been used to determine second- and third-order elastic constants. 

Energy: 

Po£(tm) = §Ciut} 3 + = f A ( v ) 

poE{r} B ) = f|Cm + C112) if + (cn + ci 2 )r? 2 = f B {rj) 
PoE(t)c) = (f Cm + IC144I rf + (5C11 + ic 44 ) n 2 = fc(v) 
PqE(t] D ) = hdn + iCies) rf + (5C11 + Icaa) = f D (v) 
PoE( VE ) = (§Ciu + 3Cu2 + C123) r? 3 + (feu + 3c 12 ) t? 2 = /b(t?) 
p E(r} F ) = C 45 67? 3 + |c 44 r? 2 = /f(j?) 

Stress: 

ti('M) = ^Cmr? 2 + cn?j = gAi(v) 
^(tja) = 5Cii2»? 2 + ci 2 r? = gA2(v) 
tslva) = (C , i23 + C , ii 2 )r; 2 +2ci27? = 5s(j7) 
Uiric) = Ci 44 r/ 2 + c 44 ?7 = gc(v) 
ts(r)D) = Ci55^ 2 + caaT) = 515(77) 

t 4 (»?F) = C 4 56?? 2 + C 44 J7 = 3f(??) 



k-point mesh) the difference between successive values 
of examined constants in our test is lower than 1 GPa. 
This difference is smaller than e.g. discrepancies observed 
between results obtained from strain-energy and strain- 
stress approach which, in the opinion of the authors, in- 
dicates that the convergence with respect to parameters 
responsible for numerical accuracy is very reasonable. 



B. Results and discussion 

Results are presented in Tableslllland lllll TableHllcon- 
tains our findings for benchmark materials Si and GaAs, 
accompanied by available experimental data and previ- 
ous theoretical findings within LDA-DFT theory. Table 
IIIII gives our prediction for the unknown values of for 
cubic nitrides. For completeness, we also provide there 
our prediction for second-order elastic moduli and com- 
pare them with previous calculations^ For c,j values, 
sometimes it was possible to determine one constant from 
a few fits [e.g., C44 from coefficients in fc(v)t fD(v)i an d 
Jf (t))]> obtaining slightly different results [e.g., for GaN, 
c 44 = 145,151,147 GPa from f c (v), fo{v), and f F ( V ) 
respectively]. In such cases the average of all obtained 
values was given in the tables. The sample plots of both 
energy and stress dependencies for GaN together with 
fitted polynomials are depicted in Figs. [5]and[3J 

When analyzing the above results, one has to bear in 
mind that both measurements and calculations of the 
third-order elastic constants are difficult. The reported 
experimental results for Cijk are determined with sig- 
nificant uncertainties and quite often exhibit discrepan- 
cies between findings of different groups (see, e.g., GaAs 
in Table [TTJ) . On the other hand, calculations of subtle 
third-order effects require reaching the limits of accuracy 
of modern quantum codes. 

When comparing experimental values with DFT re- 
sults, it is also worth noticing that ab initio calculations 
are strictly valid for perfect crystalline structure and in 
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FIG. 1: Sample convergence tests for the third-order elastic 
constants in zinc-blende GaN. Panel (a) illustrates the de- 
pendence of Gin and Gi 44 on energy cutoff (Monkhorst-Pack 
sampling llxllxll was applied for all points). Panel (b) 
shows the analogous dependence on the density of k-points 
mesh (energy cutoff 800 eV was used for all points). Note 
different scales for Cm and Ci 44 . Later on, all calculations 
for GaN have been performed with the energy cutoff 800 eV 
and the Monkhorst-Pack sampling 11 x 11 x 11. 
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FIG. 2: Plots of strain-energy relations for GaN. Squares de- 
note results of DFT computations, solid line represents a poly- 
nomial fit. See main text for details and Table|T]for definitions 
of fA{ri)J B {ri),. . . Jf{i])- 



FIG. 3: Plots of strain-stress relations for GaN. Circles denote 
results of DFT computations, solid line represents a polyno- 
mial fit. See main text for details and Table U for definitions 
of 9Ai(r)),gA2(r)) r ■ ■ , ffF (??)). Curves for g c (v) and g F (rj) co- 
incide because of very similar values of CW and C456 in GaN. 



the limit of OK temperature. The experiments, however, 
are often performed in conditions far from this idealized 
case. Particularly, the importance of temperature fac- 
tor can be verified when comparing the results of mea- 
surements for C'ijk of Si in temperatures T = 298K and 
T = 4K given in Table EE] (see Ref. [H for detailed exper- 
imental study). One can observe that for this semicon- 
ductor the values of constants C144 and C123 even change 
their sign, when the material is cooled down. 

As far as calculations of third-order elastic moduli are 
concerned, they also pose a difficult test to ab initio 
methods. The determination of is sensitive to er- 
rors in energy and stress tensor and requires extremely 
good convergence of parameters governing the accuracy 
of computations, which we believe has been reached in 
our calculations (see Fig. 1). The usage of PAW for- 
malism chosen to solve Kohn-Sham equations seems also 
not to influence the results significantly, since it has been 
demonstrated that properly performed calculations of the 
static and dynamical properties for broad range of solids 
within the PAW, pseudopotential, and LAPW schemes 
lead to essentially identical results^ In our opinion, the 
main problem lies in the approximations to the exchange- 
correlation functionals employed in various calculations. 



In the present calculations we use GGA-PBE exchange- 
correlation functional that is commonly believed to be 
one of the best in the market. However, even for the 
second-order elastic constants for GaAs (see Table UH , 
one observes significant differences between the calcu- 
lated and measured values. One possible origin of these 
discrepancies might be the commonly known tendency of 
calculations based on GGA functional to underestimate 
binding strength, and therefore to overestimate lattice 
constant. Indeed, our calculations predict the equilib- 
rium lattice constant of GaAs to be 5.75 A, considerably 
larger than the experimental value of 5.65 AM- This is op- 
posite to the local density approximation (LDA), which 
overestimates the binding and leads to lattice constants 
smaller than experimental. 

Keeping all the above in mind, we find that the agree- 
ment between our computations and measurements for 
test cases Si and GaAs is reasonably good (see Table HTl for 
details). It is also important to note that values of Cyfc 
calculated both from strain-energy and strain-stress rela- 
tions are consistent with each other. As a cross-check we 
additionally verified our approach by calculating second- 
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order elastic moduli for GaAs with the aid of the MedeA 
package, 42 It uses its own methodology of calculating Cy 
on the basis of stress computed by the VASP code.— We 
obtained values c\\ — 99 GPa, en — 41 GPa, C44 = 51 
GPa, which are in agreement with the results given in 
Table El 

Next interesting issue is to examine for which range of 
deformations the third-order effects really matter. In Fig. 
2] we compare energy and stress for the particular defor- 
mation t]b in GaN crystal with energy and stress values 
obtained within linear and nonlinear elasticity theories. 
One can clearly see that the linear approach is not suf- 
ficient for strains larger than approximately 2.5%. It is 
also worth noting that for all studied semiconductors and 
examined range of deformations (i.e., with Lagrangian 
strains up to 8%) including the terms up to third-order 
in energy expansion [Eq. |3])] sufficed to obtain good 
agreement with DFT results. 

It is also important to note that a quadratic term in e 
in the expression for Lagrangian strain rj [see Eq. (j4}] is 
usually neglected when the second-order elastic constants 
are determined. For the third-order elastic constants, 
such omission is completely unjustified. For example, the 
approximation 77 e leads to the following third-order 
elastic constants for Si, C™ ng = —256 GPa, C™ ng = 
-375 GPa, CV™ 118 = 94 GPa, C™° ng = -130 GPa, 
C^™ 118 = -105 GPa, and Cj^ ng = -2 GPa, which show 
significant disagreement with the results obtained with- 
out the aforementioned simplification (compare results in 
Table ITT]) . As one would expect, the second-order elastic 
constants remain virtually unaffected by the approxima- 
tion rj w e, now being en = 150 GPa, C12 = 62 GPa, and 
c 44 = 73 GPa. 
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FIG. 4: Energy (a) and stress component t$ (b) as a function 
of linear strain parameter e for particular deformation t\b for 
GaN cubic crystal. Full points denote results of DFT com- 
putations, solid and dashed lines indicate the curves obtained 
from nonlinear and linear elasticity theory respectively. 



IV. RELATION TO PRESSURE DEPENDENT 
ELASTIC CONSTANTS 



In the case of materials under large hydrostatic pres- 
sure it is useful to describe the nonlinear elastic prop- 
erties using the concept of pressure dependent elastic 
constants Cij(P). For many applications, it is sufficient 
to consider only terms linear in the external hydrostatic 
pressure 



(P) 
(P) 



-12 



-44 



Cll + c'uP, 
C12 + c' 12 P, 

C44 + c'44P, 



(9) 



with pressure derivatives c'ij being material parameters. 
Naturally, the information about c'ij can be recovered 
from third-order elastic constants. The necessary formu- 



las are given below-i 

2C112 



c'11 



C'l2 



Cm + 2ci2 + 2cn 



23 



2ci2 + en 
2C112 — C12 



Cll 



C 44 = 



2ci2 + cn 

2Ci55 + C144 + C44 



(10) 



2ci2 + Cn 



2C12 



Cll 



Results for pressure derivatives c'ij calculated on the ba- 
sis of our estimates for second- and third-order elastic 
constants are shown in Tables IIVI and [V] 

Table IIVI provides comparison with experimental re- 
sults for Si2£ and GaAs 4^ The agreement is very good 
and shows that the results from the strain-energy relation 
reproduce the experimental values slightly better than 
findings based on strain-stress formula. 

Table IVl contains values of c'ij for zinc-blende nitrides 
and compares the present calculation with our previous 
workfi 2 - In Ref. [l^, the following approach for the deter- 
mination of pressure dependence of the second-order elas- 
tic constants has been used. First, the hydrostatic strain 
(corresponding to the external pressure P) has been ap- 
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TABLE II: Comparison of the calculated second- and third-order elastic constants for Si and GaAs with the experimental 
values and previous calculations. All data are in GPa. 





Present results 
strain-energy strain-stress 


Previous calculations 


Experiment 




Si 




















Cll 






153 




153 


109 






lb i 


Cl2 






65 




57 


01 






00 


C44 






73 




75 


ft^ a 
OO 






sn c 

OU 


Cm 






-698 




-687 


-750 a 


-880 14 


-834 e 


-825 ' 


C112 






-451 




-439 


-480 " 


-515 d 


-531 e 


-451 ' 


C144 






74 




72 




74 


-95 c 


12' 


C155 






-253 




-252 




-385 d 


-296 e 


-310 y 


C123 






-112 




-92 





27 d 


-2 e 


-64 y 


C456 






-57 




-57 


-80 


-40 d 


-7 e 


-64 ; 


C144 + 2C155 




-430 




-432 


-580 11 


-696 d 


-687 c 


-608 ' 


GaAs 




















Cll 






100 




99 


Izo 






llo 


C12 






49 




41 


00 








l"44 






52 




51 


61 6 






60 9 


Cm 






-561 




-561 


-600 6 


-675 ft 


-622 1 


-620 > 


C112 






-337 




-318 


-401 6 


-402 h 


-387 1 


-392 3 


C144 






-14 




-16 


10 6 


-70 ,l 


2 1 


8 J 


C155 






-244 




-242 


-305 b 


-320 h 


-269 1 


-274 > 


C123 






-83 




-70 


-94 6 


_ 4 " 


-57 1 


-62 1 


C456 






-22 




-22 


-43 6 


-69 h 


-39 1 


-43 ' 


"Reference 
'Reference 
c Reference 
''Reference 
c Reference 
■'Reference 
9 Reference 
h Reference 
"Reference 
^Reference 


I (LDA). 
1(LDA). 
3$(T = 73K). 
53 (T = 4K). 
3J (T = 298K). 

II (T = 298JC). 

31 (extrapolation to T = 0-K"). 
|| ( T = 298-RT). 
.ill (T = 298K). 
13 (T = 298A"). 



plied to a crystal, and then the crystal has been addi- 
tionally deformed to determine the pressure dependent 
clastic constants. The DFT results for the total elastic 
energy combined with the strain-energy relation have en- 
abled us to determine Cij(P) as well as c\j. We would 
like to stress that the additional noninfinitesimal strain 
has not always been trace-free just leading to a spurious 
hydrostatic component that has modified external hydro- 
static pressure. Therefore, we believe that the approach 
employed in the present paper is not only more direct, 
but also slightly more accurate. The discrepancies be- 
tween our present and previous results can also be partly 
ascribed to the methodological differences, such as differ- 
ent exchange-correlation functional used and slightly dif- 
ferent calculation parameters (Brillouin zone sampling, 
energy cutoffs etc.). 



V. RELATION TO GRUNEISEN CONSTANTS 
OF LONG- WAVELENGTH ACOUSTIC MODES 

The mode Griineisen constants constitute a group of 
important coefficients, which characterize anharmonic 
properties of crystals. These quantities are frequently 
encountered in theory of phonons and in the descrip- 
tion of thcrmodynamical properties of solids. The mode 
Griineisen constants are defined as follows: 



7(qJ) 



dln^(qj) 
dlnV 



dV 



w(q,JJ 



(11) 



where u> denotes the frequency of phonon with wave vec- 
tor q and polarization vector j. V stands here for volume 
of the crystal. 

On the basis of continuum limit, one may express 
mode Griineisen constants for long-wavelength acoustic 
modes in terms of second- and third-order elastic con- 
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TABLE III: Theoretical predictions for the third-order elastic constants of zinc-blende nitrides - A1N, GaN, and InN. The 
second-order elastic constants are included and compared with previous calculations. All data are in GPa. 





strain-energy 


Present results 


strain-stress 


Previous calculations a 


a ii\r 










Cll 


Zo4 




zoz 


ZD ( 


Cl2 


ID ( 




1 AO 

14y 


1 1 1 

141 


C44 


181 
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-1070 




-1073 




Cll 2 


-1010 




-965 




C144 


63 




57 




C155 


-751 




-757 




Cl23 


-78 




-61 




C456 


-11 




-9 




GaN 










Cll 


zoo 




OKI 

ZOZ 


9^9 


C12 


1 A7 
14 / 




izy 


loi 


C44 


148 




147 


146 




-1209 




-1213 




C112 


-905 




-867 




C144 


-45 




-46 




C155 


-603 




-606 




C123 


-294 




-253 




C456 


-48 




-49 




InN 










Cll 


lOU 




loy 


14y 


C12 


110 




1 no 
lUz 


O 1 


C44 


78 




78 


77 


Cm 


-752 




-756 




C112 


-661 




-636 




C144 


16 




13 




C155 


-268 




-271 




C123 


-357 




-310 




C456 


14 




15 





"Reference Qj (GGA). 



stants. The necessary expressions used here have been 
given by Mayer and Wehner4£ The results for 7(q, j) ob- 
tained from our strain-energy estimates of elastic moduli 
are given in Table IVII 

Comparison with the experimental data available for 
Si shows that results calculated by us often differ sig- 
nificantly from experimental findings. The discrepancy 
is particularly pronounced for transverse modes (i.e., 
7 ((e,0,0),TA) = 7 ((e, e ,0),TA z ) and 7 ((e > e l 0),TA xy )) 
for which the magnitudes of Griinciscn constants are 
much smaller than for longitudinal modes. In our opin- 
ion, this indicates that 7 (q, j) are quite sensitive to inac- 
curacies in Cijk values. Therefore, one has to treat our 
prediction for mode Griineisen constants in zinc-blende 
nitrides rather as a quite crude approximation. Neverthe- 



less, it could be an interesting subject of further studies 
to compare the above results with ab initio phonon cal- 
culations performed via density functional perturbation 
theory. More detailed experimental studies for a broader 
range of materials could also shed more light on the value 
of the presented theoretical predictions. 



VI. CONCLUSIONS 

We have presented a detailed ab initio study of third- 
order elastic constants C^k for selected semiconductors 
- Si, GaAs, and zinc-blende nitrides A1N, GaN, and 
InN. Even though third-order effects are very subtle, we 
showed that it is possible to estimate them by means of 
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TABLE IV: Pressure derivatives of second-order elastic con- 
stants for Si and GaAs calculated on the basis of Eqs. (|1U[I . 
For comparison experimental findings are included. 





Present results 
strain-energy strain-stress 


Experiment 


Si 








c'n 


4.09 


4.30 


4.19 " 


C'l2 


4.34 


4.43 


4.02 


c'44 


0.27 


0.34 


0.80 


GaAs 








c'n 


4.71 


5.06 


4.63 b 


C'l2 


4.56 


4.67 


4.42 6 


c'44 


1.27 


1.48 


1.10 6 



"Reference [3| (T=4K). 
^Reference^ (T=298K). 



TABLE V: Prediction of pressure derivatives of second-order 
elastic constants for zinc-blende nitrides A1N, GaN, and InN 
calculated on the basis of Eqs. (j 10p . For comparison, results 
of previous calculations employing different methodology are 
included. 





Present results 
strain-energy strain-stress 


Previous calculations a 


A1N 








c'n 


3.53 


3.68 


5.21 


c'12 


4.12 


4.17 


4.26 


c'44 


1.03 


1.20 


1.69 


GaN 








c'n 


4.03 


4.28 


4.17 


c'12 


4.56 


4.64 


3.50 


c'44 


1.01 


1.18 


1.12 


InN 








c'n 


3.89 


4.15 


4.58 


c 12 


5.00 


5.08 


4.37 


c'44 


0.13 


0.24 


0.66 



"Reference 



for Si and GaAs with available experimental findings. 
The agreement is reasonable, however, particularly for 
moduli of smaller magnitude (e.g., for examined cases 
typically C144 and C456) relative differences are signifi- 
cant. In our opinion, they can be ascribed to three main 
factors: shortcomings of GGA-DFT theory, lack of tem- 
perature effects in our calculations (experimental results 
for Cijk are usually obtained in room temperature), and 

TABLE VI: Griineisen constants 7(q,j) for long-wavelength 
acoustic modes. Experimental results for Si were given in Ref. 
40| on the basis of ultrasound measurements data from Ref. 
291 . Theoretical prediction for 7(q,j) were calculated on the 
basis of strain-energy values for and Cijk- 





Experiment 
Si 


Theory 
Si 


A1N 


Theory 
GaN 


InN 


q= (e,0,0) 












7 (LA) 


1.108 


1.098 


1.115 


1.279 


1.415 


7 (TA) 


0.324 


0.006 


0.423 


0.459 


-0.055 


q = (e,e,o) 












7 (LA) 


1.109 


0.999 


1.066 


1.226 


1.218 


7(TA xy ) 


-0.049 


-0.301 


-0.684 


-0.613 


-1.771 


7(TA Z ) 


0.324 


0.006 


0.423 


0.459 


-0.055 


q = (e,e,e) 












7 (LA) 


1.081 


0.973 


1.056 


1.214 


1.173 



measurement uncertainties. We have also underlined the 
relation of third-order elastic constants to other anhar- 
monic properties. On the basis of the ab initio results 
for C^k, we have computed the pressure derivatives of 
second-order elastic moduli and provided rough estima- 
tions for Griineisen constants of long-wavelength acous- 
tic modes. We believe that DFT estimates of third-order 
elastic constants can be a very useful tool in modeling 
semiconducting nanostructures, in which nonlinear ef- 
fects often play an important role. 
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